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Abstract 

Results of an ongoing effort to quantify the role turbulence in scattering sound in jets are reported. 
Using a direct numerical simulation database to provide the flow data, ray paths traced through 
the mean flow are compared with those traced through the actual time evolving turbulent flow. 
Significant scattering by the turbulence is observed. The most notable effect is that upstream 
traveling waves that are trapped in the potential core by the mean flow, which acts as a wave 
guide, easily escape in the turbulent flow. A crude statistical estimate based on ray number density 
suggests that directivity is modified by the turbulence, but no rigorous treatment of non-uniformities 
in the high-frequency approximation is attempted. 


Nomenclature 

a Sound speed 

N Number of rays observed 

N r Number or rays released 

M Mach number 

p Pressure 

Re Reynolds number 

r Radial coordinate 

t 0 Jet nozzle radius 

s Entropy 

t Time 

Ui Cartesian velocities (i = 1,2,3) 

x Axial coordinate 

Xi Cartesian coordinates 

a Directivity angle 

<f) Wave phase 

p Density 

0 Cylindrical polar coordinate 

t j Angular frequency 

Subscripts 

j Jet exit 

oc Ambient 

Accents 

Q Base flow (potentially time dependent) 
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1 Introduction 

Lighthill’s theory of aerodynamic noise is oftencriticized because it does not distinguish refraction 
from generation. Instead, these effects are grouped together in a nominal ‘source’: 1 
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sound propagation ‘source’ 

where Ta is the Lighthill stress tensor. Lighthill was, of course, aware of this but concluded that 
refraction “may affect finer details, but it does not appear to be fundamental.” 2 It has since 
been argued that a distinction might not be necessary when developing predictive models, 3 but 
since generation and propagation are different physical processes, it is attractive to model them 
separately. Lighthill was not faced with the stringent noise regulations that we are today. With 
a great effort underway to achieve as little as a 3dB noise reduction, “finer details” are now more 
important. More recent but more complex acoustic analogies, such Lilley’s equation, 4,5 
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attempt to better separate propagation and generation. 

Unfortunately, since (2) is nonlinear it must be linearized for implementation and interpreta- 
tion. 5 An estimate of the steady mean flow is typically chosen to linearize about and many ongoing 
modeling endeavors take this approach. 6,7 Even modelers that circumvent the exact governing 
equations as a starting point and designate a relatively ad hoc acoustic source 8 choose to linearize 
about a steady mean flow. To linearize the propagation operator in (2), the nonlinear terms can 
be omitted, which assumes they play no substantial role at all, or moved to the right hand side, 
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Figure 2: (a) Instantaneous axial velocity: 9 evenly spaced contours from O.laoo to 0.9doo- (b) 
Instantaneous radial velocity contours: -0.3aoo to 0.3aoo with 0.1 spacing. Negative contours are 
dashed and the zero contour is omitted. 
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which is satisfying because the equation remains exact but once again blurs the distinction between 
source and propagation. Though convenient, linearization about the mean flow is well understood 
to be artificial since no individual sound wave actually encounters the mean flow. 9 Because local 
turbulence intensities can be over 100% in a jet, with large flow structures on the scale of the local 
jet radius, scattering by the turbulence might indeed be significant. We investigate this possibility. 

Refraction in jets has been investigated on many fronts. In Mach 0.5 and 0.9 jets it was studied 
by MacGregor et al . 10 by adding an artificial noise source into the jets. Suzuki & Lele 11 used 
numerical methods and analysis to study both directional and frequency scattering of sound by 
instability waves in a two-dimensional mixing layer and showed a significant influence in some 
cases. While low frequency components of jet noise typically follow the Light hill/Ffowcs Williams 
/ oc (1 — M c cos0)“ 5 law, higher frequencies are less directive, 12 fit better by three inverse Doppler 
factors. 13 High-frequency solutions 14 of Lilley’s equation for uni-directional transversely sheared 
flow might explain this, but scattering by turbulence has been offered as alternative explanation. 13 
Similar high-frequency formulations have been incorporated into predictive tools 15, 16 using a general 
high-frequency Green function derived by Durbin. 17 Adjoint Green functions have been proposed 
to simplify implementation of flow-acoustic interactions in models, 18 

The purpose of this ongoing study is to estimate the role of scattering by the unsteady turbulent 
fluctuations in a jet, and here we present a preliminary report on the effort. Since turbulence is 
analytically intractable, we rely on an existing, well-validated direct numerical simulation database 
to represent the turbulent jet flow. Flow-acoustic interaction is studied for high-frequency noise 
using an unsteady geometrical acoustic formulation to identify ray paths. Directivity is estimated 
statistically based on the ray paths, but no attempt has yet been made to construct instantaneous 
intensity profiles. 

2 Simulation Database 

Details of the direct numerical simulation database used in this study are reported in full else- 
where. 19-21 In summary, it is of a Reynolds number 3600, Mach number 0.9, temperature ratio 
Tj/Too = 0.86 turbulent jet, which matches the experimental conditions studied by Stromberg et 
al. 22 Contours of mean streamwise velocity are shown in figure 1 and contours of instantaneous 
streamwise and radial velocity are shown in 2. At this Reynolds number, the initial shear layers are 
laminar as expected and thus qualitatively different from a high- Reynolds-number jet, which would 
have turbulent shear layers, However, after transition, which occurs a little before the potential 
core closes at x % 14r 0 , the jet’s development agrees well with those at a much higher Reynolds 
number data. Downstream of the potential core, the jet’s spreading rate and Reynolds stresses 
agree with those of much higher Reynolds number jets. 21,23,24 This is important because it sug- 
gests that the energetic large scales are similar to those in jets at higher Reynolds numbers which 
should generalize the present results. 

Other points of validation are reported elsewhere. 21 In these references it was shown that the 
mean flow, the noise directivity, and far-field noise spectrum are all in excellent agreement with the 
data of Stromberg et al 22 

3 High-Frequency Approximation 

High-frequency asymptotics will be used to study the interactions of sound with the flow. This 
approach is both convenient because it leads to tractable formulations, and important because the 
high frequencies are a particularly annoying component of the noise. Of course, the noise from a jet 
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at Be = 3600 is relatively narrow banded compared to jets at typical engineering Reynolds numbers. 
However, if we accept that the large, energy carrying turbulence scales are realistic, high-frequency 
noise sources can be artificially added to the flow in order study flow-acoustic interactions. 

It is often found that high-frequency approximations give reasonable estimates for Helmholtz 
numbers as low as unity, 2 ’"* 22 and we will use this to estimate a lower bound on the Strouhal 
numbers that might be accurately represented by our procedure. We take as our length scale 
the 50 percent two-point velocity correlation width of the turbulent eddies, which, depending on 
the location within the jet, is as low as i « 0.5r o . 21 This is smaller than the scale over which 
the mean flow varies (figure 1) except in the initial shear layers. Setting the Helmholtz number 
He = ui/a = 1 gives u ~ 2a/r 0 or St = fD/Uj = Lj2r 0 /2irUj > 1.0. Of course, the motion of 
the turbulent structures might decrease or increase their effective size for a particular sound wave 
they encounter. Assuming M c ~ 0.5, this would potentially increase the Strouhal number limit by 
about one-third to 1.3, which is high but still relevant for many applications. 

We develop an unsteady ray tracing formulation similar to that used by Colonius et al. 2% We 
start with the Euler equations in Cartesian coordinates and three space dimensions, with the energy 
equation written in terms of entropy, 
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where () terms will be obtained from the simulation database and in general are functions of all 
space coordinates and time. The ()* terms are perturbations, but not necessarily acoustic at this 
point. Retaining only linear terms in the perturbations gives 
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We next assume harmonic fluctuations 


u 't 


'Ui' 

p' 


R 

p' 


P 

j'_ 


s 


gtu i<t>{x x ,t)^ 


( 7 ) 

( 8 ) 
( 9 ) 


( 10 ) 


5 


where U t , R . P and 5 are complex amplitudes that are functions of (xi,£), <f> is a real-valued phase 
function, and uj is the (large) angular frequency of the disturbances. Substituting and retaining 
only the highest order terms in uj yields 
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This leaves 5 equations and 5 unknowns. For a non-trivial solution the determinant of this system 
must be zero: 



where G = ^ 4- It is obvious that G 3 = 0 roots correspond to vorticity and entropy 

modes and are therefore not of concern to us. The remaining equation can be solved by method of 
characteristics to yield a system of ordinary differential equations for the ray paths, 
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where subscripts on <f> indicate partial differentiation. 

Once an initial coordinate is chosen, initial conditions for <f> Xi are determined by solving (15) 
with (pt = r 0 /aoo and 
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The results are insensitive to the choice of (fit and, as we shall see, (17) in general. 


4 Procedure 

Equations (15) through (18) were integrated using a fourth-order Runge-Kutta algorithm with 
A t — O.Olro/floo. Coefficients u t and a were taken from the direct numerical simulation database. 
This data was available every At = 0.17r o /aoo and at every other mesh point of the original 
computation. It was interpolated in space using B-splines and in time using a linear meth od. 
Differences were computed using second-order centered finite differences with S = 0.001r o . Results 
were insensitive to this value and did not change if a fourth-order finite difference was used instead. 
Results were also insensitive to the order of the B-spline used. We traced N r = 1000 rays through 
the mean-flow and N r — 500 rays through the unsteady flow. They were initially directed at evenly 
spaced a-angles in the 9 = 0, n and 9 = ±tt/2 planes (see figure 3). Results shown here were 
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Figure 4: Schematic showing source points used. See also table 1. 


insensitive to the number of rays traced. Except when noted, rays were traced for a total time of 
25.5r 0 / Uqo • 

Focusing by the turbulence will cause ray path to cross at points in their trajectories. At these 
caustics, the high-frequency approximation for amplitude, which can be obtained from the next 
highest terms in a/, fails because the wave fronts develop cusps, a feature that is always small with 
respect to a wavelength of the sound. In addition, behind the caustics multiple rays pass through 
each point and potentially interfere with each other. Even in a steady case with simply defined 
caustics, a rigorous solution in regions of multiple rays adds considerable complexity. In the present 
case where caustics are numerous, transitory and impossible to anticipate, we have not undertaken 
the task of developing a uniformly valid amplitude procedure. Instead, to estimate the effect of 
refraction on directivity we simply count the number of rays passing through different regions on 
a sphere surrounding the source. The resulting ray number density, N , can be shown inversely 
related to ray tube area, A, and thus proportional to tube intensity, /, which is valid where there 
axe no caustics. We assume that the volume of the caustic-affected region crossing the observation 
sphere is negligible and the positions of the caustics ever changing, as is indeed observed, and can 


Case 

x/r Q 

r/r 0 

A 

7 

0 

B 

17 

0 

C 

11 

1 

D 

17 

1 

E 

23 

1 


Table 1: Source points (see figure 4). All have 0 = 0. 
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Figure 5: Figure continues on the following page 
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(d) Case D 




Figure 5: (a-e) Bin number density for rays traced through mean and turbulent 

flows. See table 1 for source points. 




Figure 6: Every tenth ray path released in the x-y plane in the mean flow for labeled cases (see 
table 1). 


therefore neglect them in the statistical average. We also assume that behind the caustics the 
rays are decorrelated due to the stochastic action of the turbulence so we also neglect interference. 
The results of these assumptions is a statistical beam-like method which provides an estimate of 
directivity. 

Ray position statistics were accumulated as they exited a sphere of radius r = 6.0r o centered 
on the jet axis closest to the source. Monitoring a sphere at r = 10r o instead of r = 6.0r o caused 
a negligible change in our results. These data were collected in 20 bins of width Act (see figure 
3). Some rays were trapped in the jet, usually in the initial shear layers that act as a wave guide, 
and did not reach this spherical surface in the run time. These rays were not counted on the 
assumption they would have dissipated or entered the nozzle. Using a run time of l^Toja^ did not 
change the directivity estimates for a < 150°. For the unsteady case, rays were traced through 39 
separate time series that were separated in initial times by at least one large-eddy turnover time. 
The turbulence evolved as the rays progressed. 

5 Results and Conclusions 

The five point-source positions studied are labeled in figure 4 and summarized in table 1. The ray 
number-density directivities for all cases are plotted in figure 5 (a-e). In figure 5 (a), where the 
source is in the laminar potential core region (x — 7r 0 ,r = 0), we see little difference between the 
steady and unsteady ray number-density directivity. However, for rays released on the jet’s axis 
just past the end of the potential core at x = 17r 0 (figure 5 b), there is a significant difference 
between the two cases. In the mean-flow case, substantially fewer rays leave the observation sphere 
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Figure 7: Every fifth ray released in the x-y plane in the unsteady flow for the labeled cases (see 
table 1). 
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Figure 8: Mean flow ray number-density directivity for trace times of 10 r 0 ja^ and 25 r 0 ja^ 

. Case E is shown. 

at angles greater than a = 70°. The reason for this is clear in figure 6 (b): the upstream traveling 
waves are trapped by the jet's shear layers. However, in the unsteady case the rays escape and 
radiate as seen in figure 7 (b). Other source points are seen to give similar results. 

The freeing of trapped rays appears to be the principle difference between the steady and j 

unsteady cases for all the source points considered. In figure 5 (e), the upstream velocity is slow j 

enough that trapped upstream traveling waves appear as a spike in the profile near a — 180° at the j 

observation time. The trapped rays must eventually escape or enter the nozzle. Whether they will s 

eventually affect the directivity is unclear, but the rate of release is slow: the steady-flow number- 
density directivity is unchanged for ray trace times of 10r o /aoo versus 25r 0 /aoo for a < 150° (see 
figure 8). If they do eventually leave the jet, they will do so differently than the rays in a realistic j 

turbulent jet. __ _ ! 

So far we have assumed that the source itself has no inherent directivity. However, if the sources j 

are convecting quadrupoles, this is not true. In figures 9 (a-e) the rays have been given an initial j 

weighting so that they would have a three Doppler factor, W” 3 = (1 — M c cosa)“ 3 , directivity 
in absence of additional refraction. The convection Mach number was taken to be M c = 0.5, ; 

which corresponds to a convection velocity of U c = 0.6 Uj. We again see significant influence of the | 

turbulence on the directivity estimates. 

It is interesting to estimate how much of the scattering is because of the presences of the 
turbulence at any instance in time versus its change in time. To this end we have plotted a 
directivity for case B with the time derivative terms of the background flow in (18) set to zero. We § 

compare with the correct equation results in figure 10, and see that there is an effect for a > 60°, . 

but not a large one. It appears that accounting for the turbulence as a succession of steady states i 

would provide most of the scattering. I 

In summary, we have made a crude estimate of the effect of turbulence on the propagation of j 

high-frequency noise in turbulent jets and found that turbulence does affect the ray number density 
on a spherical surface outside the jet, increasing it in most cases. The most significant mechanism 
observed was that the turbulence frees rays that would otherwise be trapped traveling upstream in j 

the potential core if only refraction by the mean flow were considered. The spatial variation of the | 

turbulence was found to be more an influence on scattering than its time evolution. ; 
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Figure 10: Directivity for CASE B, x = 17 r a ,r = 0, with time derivatives in (17) and without 


This work is ongoing and a rigorous means of computing directivity as well as more realistic 
sources are being considered. 
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